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Heat transport in the cuprate superconductors YBa2Cu30 H and La2-iSr I Cu04 was measured 
at low temperatures as a function of doping. A residual linear term kq/T is observed throughout 
the superconducting region and it decreases steadily as the Mott insulator is approached from the 
overdoped regime. The low-energy quasiparticle gap extracted from kq/T is seen to scale closely 
with the pseudogap. The ubiquitous presence of nodes and the tracking of the pseudogap shows 
that the overall gap remains of the pure d-wave form throughout the phase diagram, which excludes 
the possibility of a complex component (ix) appearing at a putative quantum phase transition 
and argues against a non-superconducting origin to the pseudogap. A comparison with superfluid 
density measurements reveals that the quasiparticle effective charge is weakly dependent on doping 
and close to unity. 
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I. INTRODUCTION 

In a d-wave superconductor, the presence of nodes in 
the gap structure imposed by symmetry leads to quasi- 
particle excitations down to zero energy in the presence 
of even small amounts of disorderiS These excitations 
are delocalizcd and carry both charge and heat. The 
most striking property of this residual normal fluid is 
its universal conduction,- whereby quasiparticle trans- 
port is independent of impurity concentration. In the 
. £h ' case of heat transport, it turns out to be a direct mea- 
sure of the low-energy quasiparticle spectrum. 4 The uni- 
$_i ' versal character of heat transport was confirmed exper- 
imentally for the cuprates YBa 2 Cu 3 0,, (YBCO) 5 and 
Bi 2 Sr 2 CaCu 2 8 (Bi-2212)£ Moreover, the residual heat 
conduction measured at optimal doping! or above, in 
the overdoped regime^ is in good quantitative agreement 
with d-wave BCS theory and the quasiparticle spectrum 
either measured by angle-resolved photoemission spec- 
troscopy (ARPES) or expected from estimates based on 
the value of T c (see for instance Ref. 7). 

In this paper, we use the well-established and robust 
connection between low-temperature heat transport and 



the energy spectrum of a <i-wave superconductor to probe 
the evolution of low-energy quasiparticles and the super- 
conducting gap as a function of doping in the cuprates. 
In going from the overdoped to the underdoped regime, 
we find that the residual linear term kq/T is finite ev- 
erywhere and decreases monotonically. As a result the 
low-energy gap grows steadily, in contrast to the su- 
perconducting T c which first rises and then decreases. 
The low-energy gap in fact follows closely the normal- 
state pseudogap, measured mostly at higher energies and 
temperatures^ The ubiquitous presence of nodes and the 
tracking of the pseudogap shows that the gap remains 
of the pure (i-wave form throughout the phase diagram. 
This excludes the possibility of a complex component 
(ix) appearing at a putative quantum phase transition 
and argues against a non-supcrconducting origin to the 
pseudogap. 



II. SAMPLES 

We performed our study on two cuprate materi- 
als: the double- plane orthorhombic material YBa 2 Cu30j, 
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(YBCO) doped with oxygen in CuO chains, and the 
single-plane material I^-^Sr^CuC^ doped with Sr 
atoms (LSCO). The four samples of YBCO used in 
the study are dctwinned, flux-grown single crystals in 
the shape of platelets with typical dimensions 1.0 x 0.5 
mm and 25 /im thick. Two of them, respectively with 
y = 6.99 and y = 6.54, were grown in a BaZr03 (BZO) 
crucible , 10 which results in crystals with very high chem- 
ical purity (99.99 - 99.995%) and a high degree of crys- 
talline perfection as compared with crystals grown in 
Y 2 03-stabilised Z1-O2 (YSZ) crucibles. The sample with 
y = 6.99 was detwinned at 250°C under uniaxial stress, 
and then annealed at 350°C for 50 days, resulting in 
CuO chains with less than 0.2% oxygen vacancies, and 
hence very close to the stoichiometric composition at 
y = 7.00. 10 This level of oxygen doping is slightly above 
that for maximal T c (93.6 K), resulting in T c — 89 K. 
The sample with y = 6.54 was similarly detwinned and 
then annealed at 84° C for 2 days followed by 60° C for 5 
days. This resulted in a highly-ordered ortho-II arrange- 
ment of oxygen atoms in CuO chains, with alternating 
full and empty chains. The other two YBCO samples, 
respectively with y = 6.95 and y = 6.6, were grown in a 
YSZ crucible, and are characterized by an impurity con- 
centration typically one order of magnitude higher. The 
oxygen vacancies in the CuO chains are not ordered for 
these crystals. The y — 6.6 sample was quenched into 
an ice water bath after annealing, resulting in a higher 
level of disorder among the oxygen vacancies and thus 
a lower T c as compared to non-quenched samples with 
similar oxygen content. 

The La2- :c Sr ;c Cu04 (LSCO) samples were all grown 
in an image furnace using the travelling-solvent floating- 
zone technique and have Sr dopings of x — 0.06 (sam- 
ples A and B), 0.07, 0.09, 0.17 and 0.20. In addition, 
a non-superconducting LSCO sample with x = 0.05 
was also measured. With the exception of a; = 0.06 
B, all samples were measured as grown. This may re- 
sult in off-stoichiometric oxygen content in the samples. 
The x = 0.06 B sample was annealed in flowing argon 
overnight at 800° C in an attempt to fix the oxygen con- 
tent. The argon annealing, however, had little effect on 
our results as both x = 0.06 samples gave the same elec- 
tronic contribution to the thermal conductivity. 

In LSCO, the hole concentration per Cu in the Cu02 
planes, p, is taken to be the Sr concentration, x. In 
YBCO, however, the relation between hole concentra- 
tion and oxygen doping y is a complicated function. As 
a result, for YBCO p is determined from transition tem- 
peratures using the empirical formula™ 
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l-82.6(p-0.16) 2 



(1) 



which is a good approximation for many cuprate 
systems. — Here we use T™ ax = 93.6 K as the transition 
temperature of optimally-doped YBCO. 

The transition temperature was determined from re- 



YBCOe.99 



r/Tco=o.ooi 




T/T c 



FIG. 1: Thermal conductivity of YBa2Cu30 y and 
La2-iSr a; Cu04 vs temperature normalized at T c . Inset: the- 
oretical calculation of the effect of impurity scattering on the 
electronic thermal conductivity of cuprates (from Ref. 17). 



sistivity measurements and defined as the temperature 
at which the resistivity has fallen to zero. Note that this 
definition of T c leads to values slightly lower than those 
determined by taking the midpoint of the resistive tran- 
sition. We find however that values of T c determined 
this way correlate well with those measured by magnetic 
susceptibility and thermal conductivity. For the YBCO 
samples, the transition temperatures are 62, 44, 93.5 and 
89 K, respectively, for oxygen doping y — 6.54, 6.6, 6.95 
and 6.99. We note that the T c for the ortho-II ordered 
y = 6.54 appears anomalously high. This enhanced 
T c (consistent with the T c determined by magnetization 
measurements in similar samples^) is thought to be due 
to oxygen coordination effects, where improved oxygen 
ordering leads to a greater number of holes doped into 
the Cu02 planes as compared to a non-ordered crystal 
with the same nominal oxygen doping>i4 For the LSCO 
samples, the transition temperatures are 5.5, 8.5, 19, 16, 
34 and 33.5 K, respectively, for Sr doping x — 0.06 (sam- 
ple A, not annealed), 0.06 (sample B, annealed), 0.07, 
0.09, 0.17 and 0.20. The T c for the x = 0.09 sample 
is anomalously low, possibly as a result of Sr inhomo- 
geneity or oxygen non-stoichiometry within the crystal. 
Although different criteria may be used for determining 
the value of hole doping level in both the LSCO and 
YBCO systems, we note that small errors in our estima- 
tion of hole concentration do not noticeably affect the 
trends observed in our thermal conductivity data. 

The low temperature thermal conductivity measure- 
ments were made in a dilution refrigerator down to 40 
mK, using the standard 4- wire steady-state method with 
two Ru02 chip thermometers, calibrated in-situ against a 
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reference Ge thermometer. Currents were applied along 
the a-axis in order to probe in-plane transport and to 
avoid contributions from the CuO chains in YBCO. Ther- 
mal and electrical contact to the samples was made using 
Ag-wire and diffused Epotek H20E Ag-epoxy pads. 



III. LEVELS OF DISORDER 

It is instructive to estimate the relative amounts of 
disorder in our samples. For the YBCO samples, the in- 
clusion of impurities during growth is greatly reduced by 
using BZO over YSZ crucibles. Microwave spectroscopy 
measurements of thermally-excited quasiparticles in the 
elastic scattering limit reveal that the scattering rate in 
the superconducting state is some 12 times greater for 
optimally-doped YBCO grown in YSZ crucibles^ com- 
pared to the slightly overdoped y = 6.99 samples grown 
in BZO crucibles^ Measurements of thermal conductiv- 
ity n at high temperature, shown in Fig. 1, lead to a 
similar order-of-magnitudc difference. In the theory of 
Hirschfeld and Putikka™, the peak observed in the ther- 
mal conductivity is due to an increase in the quasiparticle 
mean free path when the sample is cooled below T c . The 
magnitude of k continues to increase with cooling until it 
becomes limited by quasiparticle scattering from impuri- 
ties and dislocations. Thus, the ratio of peak height to 
normal-state value in k(T) directly reflects the amount of 
disorder present in the crystal. The inset in Fig. 1 shows 
theoretical curves that demonstrate this effect fli where 
the electronic contribution to thermal conductivity nor- 
malized by the value of k at T c is plotted as a function 
of T/T c . " 

Samples with a large ratio of impurity scattering rate 
(given by T) to T c q are seen to have a peak height that 
is suppressed. An order of magnitude increase in the in- 
trinsic level of disorder within the crystal results roughly 
in a factor of two decrease in the peak height. The data 
presented for our crystals in Fig. 1 reflects that effect. 
The deliberate addition of impurities such as Zn in YSZ- 
grown optimally-doped YBCO samples leads to a large 
suppression in peak height. Data for a sample with 1.7% 
Zn impurities (determined from T c suppression) is shown 
in Fig. 1, and it is seen that the addition of this level of 
impurities causes the peak to nearly vanish. The corre- 
sponding residual resistivity extrapolated from the linear 
temperature dependence of p a (T) goes from being nega- 
tive in the nominally pure crystal to po = 30 pQ cm in 
the Zn-doped crystal. It is clear that the optimally-doped 
LSCO sample (x = 0.17) shown in Fig. 1, with po = 33 
pfl cm (as extrapolated from resistivity data above T c ), 
exhibits much stronger scattering than any of the YBCO 
samples. This is true despite the high chemical purity 
of the crystal, and is likely a result of the Sr atoms in- 
cluded as dopants acting also as scatterers. Considering 
all the evidence we estimate the relative amount of disor- 
der in the various crystals studied here to be roughly in 
the proportion of 100:10:1 for LSCO, YSZ-grown YBCO 
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33.5 


0.20 


330 ± 40 


36 




Bi-2212 


89 


0.16 


150 ± 30 


19 


30 


Tl-2201 


15 


0.26 


1400 ± 70 


270 
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TABLE I: Compilation of T c , doping and residual linear term 
in the thermal conductivity as well as values of the quasipar- 
ticle anisotropy ratio vf/v2 from Eq. 2 and gap maximum 
Ao (see caption of Fig. 6) for the samples in this study. Data 
for optimally doped Bi-2212^ and overdoped T1-2201& from 
previous studies are provided for completeness. 

and BZO-grown YBCO, respectively. 



IV. DOPING DEPENDENCE OF k /T 

The low-temperature thermal conductivity of YBCO 
and LSCO samples is shown as a function of tempera- 
ture in Figs. 2 and 3, respectively. The data is plotted 
as k/T vs T 2 because the quantity of interest is the resid- 
ual linear term k /T, defined as the T — limit of k/T, 
obtained by extrapolation of the low-temperature data. 
This residual linear term can only be due to fermionic 
carriers and is attributed to zero-energy quasiparticles. 
Indeed, as will be seen below, it is a direct confirmation, 
via a robust bulk measurement, of the d-wave nature of 
the superconducting order parameter in cuprates. The 
extrapolation procedure is described in detail in the Ap- 
pendix, where the contribution of phonons is analyzed. 
The main results of the paper do not depend on the 
particular extrapolation procedure. This is true, for ex- 
ample, for the overall trend with doping, which is im- 
mediately evident from Figs. 2 and 3: k /T decreases 
steadily with underdoping, all the way from the slightly 
overdoped to the highly underdoped regime. Using the 
extrapolation procedure outlined in the Appendix, the 
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FIG. 2: Thermal conductivity of YBa2Cu30 y plotted as k/T 
vs T 2 . Open symbols represent samples grown in YSZ cru- 
cibles and filled symbols those grown in BZO crucibles. 



values we obtain are given in Table I. Note that a mea- 
surement on a fully deoxygenated YBCO sample with 
y = 6.0 correctly yields a zero intercept: k /T = ± 1 
/iW K -2 cm -1 . The values for LSCO agree with those 
published in a previous study«iS, with the exception of 
their x = 0.17 sample which has been measured to be ap- 
proximately twice the value we observe. We attribute this 
difference to the fact that the crystal studied by Takeya 
et al. had a T c of 40.2 K compared to our T c of 34.2 
K, pointing to a slightly higher hole concentration (likely 
due to different oxygen levels within the crystals). 

Let us analyze these results within the framework of 
standard d-wave BCS theory. In the clean limit at low 
temperature, when fcgT <C 7 <C ksT c , where 7 is the im- 
purity bandwidth, the quasiparticle thermal conductivity 
can be written as£ 
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where n is the number of Cu02 planes per unit cell and 
d is the c-axis lattice constant, vf and vi are the quasi- 
particle velocities normal and tangential to the Fermi 
surface at the node, respectively, and are the only two 
parameters that enter the low-energy spectrum, given by 
E = h^Jvpk\ + v\k\ where ki and k 2 are vectors 
normal and tangential to the Fermi surface at the node, 
respectively. The parameter vi is simply the slope of the 
gap at the node: 
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FIG. 3: Thermal conductivity of La2- ;E Sr; C Cu04 plotted as 



k/T vs T\ for a) 



0.09, 0.17 and 0.20, and b) x = 0.05, 



0.06 and 0.07. The lines through the data are power law fits, 
discussed in the Appendix. 



where kp is the Fermi wavevector at the nodal position. 
These are remarkably simple formulae, which provide a 
direct access to the parameters that govern low-energy 
phenomena in a e?-wave superconductor. The residual 
heat conduction in Eq. 2 is not only universal, i.e. in- 
dependent of scattering rate (or impurity bandwidth), 
but it was also shown to be independent of Fermi-liquid 
corrections and vertex corrections (i.e. corrections due 
to anisotropic scattering between nodes) A Note, how- 
ever, that the latter two corrections affect the microwave 
(charge) conductivity (see below). 

In Fig. 4, the anisotropy ratio vpjvi is plotted against 
carrier concentration p, using Eq. 2 and the values of 
kq/T listed in Table I. Also included is the published 
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FIG. 4: Anisotropy ratio Vf/vz, calculated from thermal 
conductivity data via Eq. 2, vs hole doping per Cu02 plane, 
p. The data for Bi-22f 2 is from Ref. 7. Lines are guides to 
the eye (solid for YBCO, dashed for LSCO) 



value for Bi-2212 at optimal dopingii All three cuprates 
have a comparable anisotropy ratio at optimal doping: 
v F /v 2 = 10, 12 and 19, for LSCO, YBCO, and Bi-2212, 
respectively. It has already been noted 7 - that the value 
of 19 for Bi-2212 is in excellent agreement with the ratio 
of 20 coming from values of vp = 2.5 x 10 7 cm/s and 
v 2 = 1.25 x 10 6 cm/s obtained directly from ARPES-. 
(Note the value of 12 for the optimally-doped YBCO 
crystal differs slightly - albeit within error bars - from 
our previously published result of 14, which was an aver- 
age of several samples £) 



V. DISCUSSION 

A. Nature of the superconducting order parameter 

Several authors have proposed the existence of a quan- 
tum critical point within the superconducting dome in 
the phase diagram of cuprates, either as a theoretical 
prediction to explain the diagram itself or as suggested 
in various experiments. Its location is usually taken to 
be near (or slightly above) optimal doping, in the neigh- 
bourhood of p = 0.2. If it is associated with a change in 
the symmetry of the superconducting order parameter, 
Vojta et al. have argued that the most likely scenario is 
a transition from a pure d x i_ y i state to a complex order 
parameter of the form d x 2_ y 2 + ix, where x can have ei- 
ther s or d xy symmetry»22, Sharoni et al. have recently 
reported a split zero-bias anomaly in their tunneling on 
Y-123 thin films as soon as the material is doped be- 
yond optimal doping, a feature which they attribute to 



the appearance of a complex component to the order pa- 
rameter in the bulk^i The presence of a subdominant 
component ix in the order parameter causes the nodes 
to be removed, as the gap can no longer go to zero in 
any direction. Our observation of a residual linear term 
in the thermal conductivity of both YBCO and LSCO, 
as well as previous results on optimally-doped Bi-2212 7 
and strongly-overdoped Tl-2201— , is a direct consequence 
of nodes in the gap. It therefore excludes the possibil- 
ity of any such subdominant order parameter in the bulk 
throughout the doping phase diagram. In other words, 
if there truly is a quantum critical point inside the su- 
perconducting dome, it does not appear to be associated 
with the onset of a complex component in the order pa- 
rameter. 

In view of the ubiquitous nature of the residual lin- 
ear term in superconducting cuprates, observed in four 
different hole-doped materials from strongly-overdoped 
Tl-2201 to strongly-underdoped LSCO, two previous re- 
sults stand out as puzzling anomalies: the absence of a 
detectable linear term in electron-doped Pr2_ K Ce I Cu04 
(PCCO) 22 and in hole-doped YBa 2 Cu 4 8 ^2, In particu- 
lar, note that the upper bound of 0.02 mW K -2 cm -1 
quoted for kq/T in YBa2Cu40s is 4 to 5 times lower than 
the value obtained here for YBa2Cu30 y at a comparable 
hole concentration (y = 6.54 or 6.6) - as assessed by the 
very similar resistivity curves above T c - and compara- 
ble sample quality. This extremely low value is akin to 
that found in non-superconducting strongly-underdoped 
LSCO (x = 0.05). 



B. Effects of disorder 

One of the most remarkable results of transport theory 
in d-wave superconductors is the universal nature of heat 
conduction, which appears due to a cancellation between 
the increase in normal fluid density and the decrease in 
mean free path observed as the concentration of impuri- 
ties is increased^ This universal behaviour is only found 
in the clean limit where HT <C Ao. In situations where T 
is large, (or Ao is small), the behavior is no longer uni- 
versal, and the measured linear term may be closer to the 
normal state value kn/T than the universal limit. In the 
extreme case where HT ~ Ao, superconductivity is de- 
stroyed and the normal state value of kn/T is recovered. 
Therefore the validity of using Eq. 2 to extract values of 
vp jv% from measurements of the residual linear term is 
ensured only when samples are in the clean (universal) 
limit, HT <C Ao. Universal behaviour in YBCO at opti- 
mal doping is already well established, 5 " and inspection of 
Fig. 4 shows that this is confirmed at other dopings. In- 
deed, we observe that both BZO and YSZ grown crystals 
yield values of Vp/v2 that lie on the same curve despite 
having an order of magnitude difference in purity level, 
which is strong evidence that the clean limit is reached 
in our YBCO samples. 

In LSCO, the extremely small values of kq/T measured 
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in highly underdoped samples point to a different conclu- 
sion. Indeed, for x = 0.06, k /T ~ 12 fiW K~ 2 cm" 1 , 
while the minimum value for LSCO allowed by Eq. 2 is 
|f f(l + 1) = 18.3 K- 2 cm" 1 . The data for the 
LSCO samples with the lowest dopings are plotted in 
Fig. 5, which shows that the use of Eq. 2 for these samples 
is invalid. This breakdown suggests that our underdoped 
LSCO samples are not in the clean limit, and hence we 
cannot extract quantitative information by using Eq. 2, 
as we will do for YBCO in the following sections. The 
same conclusion would apply to previous LSCO data—. 

In order to understand the LSCO data within a d- 
wave BCS theory of low temperature heat transport, it 
will be necessary to incorporate the effects of impurity 
scattering in the underdoped regime outside of the clean 
(universal) limit. The effect of impurity scattering on a 
d-wave superconductor has been worked out in the stan- 
dard case of a normal state that is metallic, and conducts 
heat better than the superconducting stated When the 
concentration of impurities is increased in such a case, 
T c is gradually suppressed to zero and the residual linear 
term rises monotonically to meet its normal state value. 
However, our LSCO samples with x < 0.09 exhibit the 
well-known insulating upturns in the normal state resis- 
tivity associated with the ground state metal-insulator 
transition observed near x ~ 0.16m^ In fact the resis- 
tivity in a strong magnetic field appears to diverge as 
T — > Oi^Si Thus, for the LSCO samples where x < 0.16, 
the effect of increasing the impurity concentration would 
be to evolve the system towards an insulating state, or at 
least one that conducts heat less well. In this scenario, 
we expect the measured residual linear term kq/T to be 
smaller than the universal value, which would explain 
how in Fig. 5 we measure a linear term smaller than that 
allowed by Eq. 2. 

Another possibility is suggested by the theoretical work 
of Atkinson and HirschfeldSI, in which the Bogoliubov- 
deGennes equations are used to model the paired state 
as an inhomogenous superfluid. This approach allows for 
the possibility of quantum interference processes such as 
localization which are neglected in the usual framework. 
In their model, the residual linear term Ko/T is seen to 
decrease in the presence of increasing impurity concentra- 
tion, a direct result of weak localization of carriers. The 
fact that we measure a linear term in underdoped LSCO 
which is smaller than that allowed by Eq. 2 may be evi- 
dence for the existence of such localization in LSCO. We 
hope these observations will stimulate further theoretical 
work. 



C. Doping dependence of the superconducting gap 

The remarkable success of Eq. 2 at optimal doping 
validates the extension of our study across the doping 
phase diagram, at least for our YBCO samples, where 
the clean (universal) limit is established. In interpret- 
ing our measurements of the anisotropy ratio Vf/v2 in 
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gap of the form Abcs = 2.14fcsT c is also plotted, with T c 
taken from Eq. 1 (and T™ ax = 90 K). The value of the pseu- 
dogap in Bi-2212, as measured by various techniques;- is also 
shown (open symbols). The thick dashed lines is a guide to 
the eye. 
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such a study, the first thing to emphasize is the fact that 
vp, the Fermi velocity at the node, is essentially inde- 
pendent of doping. This was shown by ARPES both in 
Bi-2213i2 and in LSCO^ where the slope of the E vs 
k dispersion at the Fermi energy is seen to vary by no 
more than 10% over the range 0.03 < x < 0.3, with an 
average value of vf ~ 2.5 x 10 7 cm/s in both materials. 
The position of the node in fc-space is also independent 
of doping^, with Uf — 0.7 A ~ 1 as measured from (tt,tt) 
to the Fermi surface. As a result, a study of k /T vs p 
yields the doping dependence of v 2 = v 2 (p). In Fig. 6, 
we plot the slope of the gap at the node as a function 
of carrier concentration, not as v 2 vs p but in a more fa- 
miliar guise as the corresponding gap maximum, Ao, of 
a putative <i-wave gap function A = Aocos20, via Eqs. 2 
and 3. Given that kp is constant, this is equivalent to 
plotting v 2 directly. The values of Ao are also listed in 
Table I. Again, here we have confined our analysis to 
YBCO only, given that LSCO was seen to lie outside the 
clean limit. Plotted alongside this data is a conventional 
BCS d-wave gap (dashed curve), where we have assumed 
Ao = 2.14fcsT c (weak-coupling approximation). The p 
dependence of the gap is estimated using Eq. 1, with a 
maximum T c at optimal doping of 90 K. 

Let us examine the implications of these results by 
starting on the overdoped side of the phase diagram. The 
only available data in the strongly-overdoped regime is on 
Tl-2201, 8 a single-plane cuprate with optimal T c ~ 90 K. 
For an overdoped crystal with T c = 15 K, the measured 
residual linear term is kq/T = 1.4 mW K~ 2 cm" 1 , which 
yields vf/v 2 = 270 via Eq. 2. In comparison, the weak- 
coupling BCS prediction based on the value of T c — 15 K 
is vf/v2 = 210, using the values of Dp and kp given 
above. The good quantitative agreement shows that in 
this strongly-overdoped regime BCS theory works quite 
well, and the much larger anisotropy ratio is a conse- 
quence of the much smaller T c . 

We now turn our attention to the undcrdoped region 
of the phase diagram. In the case of YBCO the decrease 
in kq/T by a factor 2 between y = 6.99 and y — 6.54 
provides one of the main results of this paper: the ve- 
locity ratio decreases with underdoping; it drops from 
16 to 8 in going from a sample with T c = 89 K to an 
underdoped sample with T c = 62 K. This reflects an 
underlying steepening of the gap at the node while T c 
drops, with underdoping. (Note that this is in contra- 
diction with the results of Mesot et al. who extracted a 
slope of the gap from their ARPES measurements on Bi- 
2212 near optimal doping that seemed to decrease slightly 
with underdoping^ and with the analysis of Panagopou- 
los et al. who extract a gap maximum from their pene- 
tration depth measurements that remains approximately 
constant in the underdoped regime^) 

Taken by itself, this could be attributed either to a 
gradual departure from weak-coupling towards strong- 
coupling BCS superconductivity, with a growing ratio 
A /T c . It could also be interpreted as a gradual defor- 
mation of the gap shape, from a simple cos2^ angular 



dependence to a much steeper function with a decreas- 
ing average gap that scales with T c . However, in view 
of the known behaviour of the pseudogap, these expla- 
nations are unlikely to be the main story. Indeed, the 
growth of the low-energy gap observed through kq /T is 
highly reminiscent of the similar trend observed in the 
high energy gap ( or pseudogap ) with underdoping. 
In fact the growth of Ao derived from v 2 is in quanti- 
tative agreement with the pseudogap maximum deter- 
mined by ARPES; 30 i 31 i 32 tunneling 3 ^ 3 ! a - b plane opti- 
cal conductivity 3 ^ and Raman scattering) 3 ^ 7 , as shown 
in Fig. 6. 

This striking similarity in scaling points to a common 
origin, which allows us to say the following things on the 
nature of the pseudogap. First, due to the very existence 
of a residual linear term, the (total) gap seen in thermal 
conductivity at T — » is one that must have nodes. Sec- 
ondly, it has a linear dispersion as in a d-wave gap (i.e. it 
has a Dirac-like spectrum). Thirdly, it is a quasiparticle 
gap and not just a spin gap. A fundamental question is 
whether the pseudogap is related to or independent of 
superconductivity. The first and most natural possibil- 
ity is that it is due to some form of precursor pairing. A 
second possibility is that it may come from a distinct non- 
superconducting state. Indeed, a universal thermal con- 
ductivity is also possible in a non-superconducting state 
as long as the energy spectrum is Dirac-like (i.e. linear 
dispersion). For example, a universal (charge) conductiv- 
ity was derived for a degenerate semiconductor in 2Dm& 
Interestingly, the d-density-wave (DDW) state proposed 
as an explanation for the pseudogap phenomena seen in 
underdoped cuprates^S, also exhibits a universal conduc- 
tivity provided that the chemical potential /i = 0. In 
the region where both orders coexist - DDW and <i-wave 
superconductivity (DSC) - Eq. 2 is then predicted to 

hold A with v 2 replaced by ^J(v% DW ) 2 + (v£ c ) 2 , where 
v ddw anc j V SC are i-j^g ve i oc ities for the two types of 
order, respectively. The main question then is how does 
the chemical potential evolove as a function of doping. 

In summary, our measurements of k/T throughout the 
phase diagram allow us to make the following statements 
about the evolution of A with doping. First, the extrap- 
olated value of the gap maximum from thermal conduc- 
tivity in the overdoped regime is in excellent quantitative 
agreement with that expected from BCS theory. Sec- 
ondly, Ao continues to grow with underdoping while T c 
rises and then falls, in contradiction to what one would 
expect from BCS theory. The divergence of these two 
energy scales in the underdoped regime is a manifesta- 
tion of the pseudogap, whose presence is now revealed 
at very low energies in a bulk measurement on crystals 
of the utmost quality and purity. The fact that the gap 
preserves its pure d-wave form (with nodes on the Fermi 
surface) throughout strongly suggests that the pseudogap 
is superconducting in origin. 
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D. Superfluid density and microwave conductivity 

One way to shed further light on the nature of the low- 
energy electron state in underdoped YBCO is to com- 
pare heat transport and charge dynamics. For a e?-wave 
BCS superconductor, Durst and Lee have shown that the 
two conductivities are affected differently by scattering 
anisotropy and quasiparticle interactions^ The charge 
conductivity in the u> — > and T — > limit is given 
by 



r (rr\ e2 1 n a 2 V F 

hm T ^o cri(T) = (T = t-o-^ Pvc a FL 



(4) 

V2 



where e is the electron charge. The factor (3yc is due to 
vertex corrections and is greater than f .0 when impurity 
scattering is anisotropic. This simply reflects the fact 
that intra-node scattering will degrade a charge current 
less than inter-node (opposite- or side-node) scattering 
that involves a larger change in momentum. This is the 
discrete version of the (1 — cosO) term that enters nor- 
mal state conductivity and reflects the predominance of 
back-scattering over small-angle scattering. Numerical 
calculations suggest that five can be large (e.g. in ex- 
cess of fO) in high-purity samples as long as impurity 
scattering is in the unitary limit A Note that vertex cor- 
rections have negligible effect on heat transport. The fac- 
tor a 2 FL is a Fermi-liquid correction which arises because 
of quasiparticle-quasiparticle interactions. The same fac- 
tor also enters in the low-temperature slope of the normal 
fluid density p n (T) = Ps (T = 0) - p.(T)A 



p n (T) _ 21n2 f 



h 2 d 



l FL 



V-2 



k B T. 



(5) 



The temperature dependence of the a-axis superfluid 
density of YBCO crystals very similar to ours was mea- 
sured via the penetration depthi 16 ! 41 The value of a 2 FL ^ 
obtained from this data is shown in Fig. 7. Using the 
value of vp/v2 from kq/T (averaging the YSZ-grown and 
BZO-grown data) yields 




FIG. 7: Quasiparticle velocity ratio in YBCO obtained from 
universal heat transport, as v F /v2 (circles) via Eq. 2, and 
from superfluid density data of Refs. 16 and 42, as a% L ^ 
(triangles) via Eq. 5. Lines are guides to the eye. 



v F , along with a V2 that increases with decreasing doping 
provides a significant challenge to microscopic theories of 
d-wave superconductivity in cuprates. 

The microwave conductivity ai(uj, T) of YBCO was 
recently measured in crystals nominally identical to ours 
with y = 6.50^i and y — 6. 99. 16 Even though the mea- 
surements go down to 1 GHz and 1.3 K, it turns out to be 
unclear how to reliably extrapolate this data to the uj = 
and T — limit, so that a meaningful comparison of kq/T 
and (Jo is not quite possible at this stage. The shape and 
temperature dependence of the microwave spectrum for 
the y = 6.50 sample for example is suggestive of non- 
unitary scattering close to the Born limit, implying that 
the low-temperature universal limit regime may not be 
reached by 1.3 K. Further work is needed to ascertain 
whether this is indeed the correct scenario. 



a 2 FL ~ 0.4 - 0.5 , at p ~ 0.16 (6) 



a 2 FL ~ 0.6 - 0.7 , atp- 0.09 (7) 

A similar value was previously derived for optimally- 
doped Bi-2212. 7 We conclude that this FL parameter is 
near unity and, more importantly is only weakly depen- 
dent on doping. In a recent paper, loffe and Millis 42 
argue that a doping independent a FL , interpreted as ef- 
fective charge, is inconsistent with both the Brinkman- 
Rice mean field theory and slave boson gauge theory ap- 
proaches to the Mott physics of high-T c materials. In- 
deed the combination of a doping independent ot 2 FL and 



VI. CONCLUSIONS 

We have studied the evolution of thermal transport as 
T — » in the cuprate superconductors YBCO and LSCO 
over a wide range of the doping phase diagram. The 
residual linear term, ^ , is observed to be finite through- 
out the superconducting region. This proves that the gap 
always has nodes on the Fermi surface, a fact that has 
two implications: 1) it rules out the possibility of a mul- 
ticomponent order parameter of the type d + ix in the 
bulk, appearing at a putative quantum phase transition, 
and 2) it argues strongly in favour of a superconducting 
origin to the pseudogap (e.g. precursor pairing). As the 
Mott insulator is approached, ^ is observed to decrease, 
leading to a decreasing value of the quasiparticle velocity 



9 



anisotropy ratio, This result offers some of the 

first insights into the doping dependence of several im- 
portant quasiparticle parameters. First, the slope of the 
d-wave superconducting gap at the nodes, v 2 , is seen to 
increase steadily as doping is decreased, consistent with 
a growth of the gap in the underdoped regime. This is in 
contradiction to what one naively expects from BCS the- 
ory, where the gap scales with T c . The gap we extract at 
very low energies follows closely the pseudogap measured 
mostly at much higher energies by other techniques. This 
close tracking of the pseudogap shows that the gap re- 
mains roughly of the pure d-wavc form throughout the 
phase diagram. Secondly, a comparison with supcrfluid 
density reveals that the quasiparticle effective charge is 
weakly dependent on doping and close to unity. 

The considerable difference between the magnitude of 
the change in ^ with underdoping for the LSCO and 
YBCO samples provides clues as to the role of disor- 
der in the underdoped regime. In particular, the small 
value of the residual linear term of the most highly under- 



doped LSCO samples is incompatible with the standard 
theory of transport for d-wave superconductors, motivat- 
ing theoretical work which would incorporate the effects 
of impurities in a superconductor whose normal state is 
insulating rather than metallic in nature. 
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FIG. 8: Thermal conductivity of the s-wave superconductor 
V3SL The data is plotted as ^ vs T 1 ' 74 , and the line repre- 
sents a free fit to the data of the form of Eq. 9. The resulting 
linear term is zero: Oil fiW K -2 cm -1 , consistent with that 
expected for a nodeless superconductor. 



VIII. APPENDIX: PHONON THERMAL 
CONDUCTIVITY IN D-WAVE 
SUPERCONDUCTORS 

In order to use thermal conductivity as a direct probe 
of low-energy quasiparticles in <i-wave superconductors, 
the contribution from phonons must be reliably ex- 
tracted. This may be achieved by performing experi- 
ments in the regime T — > 0, where the phonon mean free 
path becomes limited only by the physical dimensions of 
the sample. From simple kinetic theory, the conductivity 
of phonons in this boundary-limited scattering regime is 
given by 

Kph = ^0 <v ph > A T 3 (8) 

where j3 is the coefficient of phonon specific heat, Ao is 
the temperature-independent mean free path, and <v p h> 
is a suitable average of the acoustic sound velocities. The 
electronic linear term is then naturally extracted by plot- 
ting thermal conductivity data as y vs T 2 and interpret- 
ing the intercept as the residual linear term at T = 0, and 
the slope as the phonon contribution governed by Eq. 8. 
The extension of our measurements into the highly un- 
derdoped region of the cuprate phase diagram, where y 
becomes very small, led us to refine this extrapolation 
technique. 

To motivate why this may be necessary, consider the 
possible scattering mechanisms available to a phonon im- 
pinging upon the surface of a crystal. The phonon may 
either be absorbed and reemitted with an energy distri- 
bution given by the local temperature (diffuse scattering) 



or it may be reflected elastically (specular reflection) . In 
the case of diffuse scattering, the phonon is reradiated in 
a random direction resulting in a temperature indepen- 
dent value of Ao and a T 3 dependence of K p h as recognised 
by Casimir^ 3 " However, as the temperature of a crystal 
is reduced and the average phonon wavelength increases, 
a surface of a given roughness appears smoother, which 
may increase the occurance of specular reflection and re- 
sult in a mean free path which varies as some power of 
temperature, so that n p h oc T a . We would thus expect 
a deviation from the diffuse scattering limit of T 3 tem- 
perature dependence of n p h for samples with sufficiently 
smooth surfaces. Such an effect has been previously ob- 
served in many studies of low-temperature phonon heat 
transport in high-quality crystals, such as AhOs^ Si^ 
KC1 and KBr4& LiFii and diamond^ 

This effect can be seen in most of our samples, mani- 
festing itself as a gradual curvature in the low tempera- 
ture part of our data when plotted as ^ vs T 2 (see Figs. 2 
and 3). In light of this, the thermal conductivity in the 
boundary scattering regime is more correctly modelled 
as: 

k = K e i i Kph = AT + BT a (9) 

with a < 3. Here A is the coefficient of the electronic lin- 
ear term, and B the temperature-independent coefficient 
of the phonon term, where a is some power of tempera- 
ture, typically between 2 and 3. (Note that there is no 
fundamental reason for a single power law - it is sim- 
ply an empirical result. For example, in AI2O3 previous 
studies™ have found a = 2.77.) 

In superconductors possessing an isotropic or s-wave 




0.0 0.2 0.4 

jl.71 ( K 1.71 ) 

FIG. 9: Thermal conductivity of YBCOe.e- The data is 
plotted as ^ vs T 1 ' 71 , and the line is a linear fit. Inset: zoom 
at low temperatures. Note the presence of a clear residual 
linear term, the contribution of nodal quasiparticles. 
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gap, the absence of an electronic linear term at low tem- 
peratures reveals this effect well. Plotted in Fig. 8 is 
thermal conductivity data for the s-wave superconduc- 
tor V^Si^ where the line is the result of a free fit to a 
simple power law as in Eq. 9. Such a procedure yields 
a linear term A = -0.04 ± 1 /xW K -2 cm" 1 , a phonon 
coefficient B = 5.73 ± 0.07 mW R-^+^cm" 1 , and an 
exponent a — 2.74 ± 0.01. The validity of such a fitting 
procedure is best seen by plotting the data as in Fig. 8, 
with the cc-axis in units of T a_1 . The striking linearity 
of the data on this plot, and the fact that it extrapolates 
to zero, is good evidence for the appropriateness of the 
power law fitting procedure. 

Fig. 9 shows the results of a fit to Eq. 9 for our 
YBC06.6 crystal, where a power law of a = 2.71 is seen 
to persist to temperatures as high as 550 mK. Power- 
law fits are also shown in Fig. 3b, this time on a k/T vs 
T 2 plot, for underdoped LSCO samples. The value of a 
observed in our samples was found to vary over a range 
from 2.4 for the YBCO y = 6.99 crystal to 2.92 for the 
LSCO x = 0.09 crystal. 

It is worth stressing that the single-power-law fitting 
procedure described here is simply an empirical approach 
to extrapolate the most reliable value of k/T at T = 0. 
As a three-parameter free fit to the data over a tem- 
perature range typically of a decade (50 - 500 mK), it 
is far better than the old two-parameter forced fit to a 
k/T = a + bT 2 form, which invariably must be limited to 
the very lowest temperatures (usually below 150 mK or 
so) and typically overestimates the value of Ko/T. How- 
ever, it must be noted that in some cases it doesn't work 



well over the whole range up to ^500 mK. This is indeed 
the case in our LSCO samples x = 0.17 and x = 0.20, 
where the single-power law fit is inadequate to describe 
the rapid fall of k/T below 150 mK. (Such a decrease was 
also observed in LSCO samples with similar doping lev- 
els in a previous studyi&.) The low temperature drop is 
most likely in the electronic channel (K e (T)), but its ori- 
gin is as yet unclear. A downturn in k/T at temperatures 
below 0.2 K or so can be induced in a sample on purpose 
by simply using highly resistive contacts (kTl or higher). 
The drop is then attributed to the rapid deterioration of 
the coupling between electrons and phonons at those very 
low temperatures. In such a case, the extrapolation pro- 
cedure must be based only on data above the downturn. 
It is not clear that the same phenomenon can still occur 
in the presence of excellent contacts, like those used here 
(less than 1 f2). These considerations are explored and 
discussed more fully elsewhere^. In conclusion, when a 
single-power law works over a wide range of temperature 
(e.g. up to 0.5 K), then the extrapolation is reliable; if 
it doesn't work, then one needs to understand why and 
may be forced to rely only on data above any anomalous 
downturns. 

In this paper, all data was successfully analyzed using 
the power-law procedure, except for LSCO samples x = 
0.17 and x = 0.20, where instead a linear fit to the data 
in Fig. 3a was used below 150 mK, yielding the values 
of Ko/T quoted in Table I. Use of a power-law fit above 
150 mK yields higher values, namely kq/T = 0.2 and 
0.4 mW K -2 cm -1 respectively, which has no impact on 
any of our conclusions. 



